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The strong coupling constant is a fundamental parameter of nature. It can be extracted from ex- 
periments measuring three-jet events in electron-positron annihilation. For this extraction precise 
theoretical calculations for jet rates and event shapes are needed. In this talk I will discuss the 
NNLO calculation for these observables. 
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1. The calculation 

The process e^e^ — )• 3 jets is of paiticular interest for the measurement of the strong coupling 
a,. Three-jet events are well suited for this task because the leading term in a perturbative calcula- 
tion of three-jet observables is already proportional to the strong coupling. For a precise extraction 
of the strong coupling one needs in addition to a precise measurement of three-jet observables in 
the experiment a precise prediction for this process from theory. This implies the calculation of 
higher order corrections. The process e+e^ — )• 3 jets has been been calculated recently at next-to- 
next-to-leading order (NNLO) in QCD ^. The master formula for the calculation of a three-jet 
observable at an electron-positron collider is 



°^ n>3 
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(1.1) 



where q\ and q2 are the momenta of the initial-state particles and 1/(8^) corresponds to the flux 
factor and the average over the spins of the initial state particles. The observable has to be infrared 
safe, in particular this implies that in single and double unresolved limits we must have 
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J's , . . . , P3 , <7i , ^2 ) for single unresolved Umits , 
h{p\,---,P3,qi,q2) for double unresolved limits . 



(1.2) 



is the amplitude with n final-state partons. At NNLO we need the following perturbative ex- 
pansions of the amplitudes: 
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Here £/n '^ denotes an amplitude with n final-state partons and / loops. We can rewrite symbolically 



the LO, NLO and NNLO contribution as 
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The computation of the NNLO correction for the process e+e^ — > 3 jets requires the knowledge of 
the amplitudes for the three-parton final state e^e^ — )• qqg up to two-loops the amplitudes 

of the four-parton final states e^e^ qqgg and e^e^ — > qqq'q' up to one-loop [§-|8|] and the five- 
parton final states e+e^ — )• qqggg and e^e^ — > qqq'q' g at tree level [^-|ll|]. The most complicated 
amplitude is of course the two-loop amplitude. For the calculation of the two-loop amplitude 



special integration techniques have been invented [12 - 15 1. The analytic result can be expressed in 
terms of multiple polylogarithms, which in turn requires routines for the numerical evaluation of 
these functions |[l|, |l7|]. 
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2. Subtraction and slicing 

Is is well known that the individual pieces in the NLO and in the NNLO contribution of 



eq. (1.4) are infrared divergent. To render them finite, a mixture of subtraction and slicing is 



employed. The NNLO contribution is written as [18] 



NNLO _ I ^^(0) _ ^ ^^^ngle _ ^ ^^(0,2) 

+ / ( ^4 do^^'^ + ^4 oc/af^'" - e^oda^^'^'^ 
+ / ( ^3 dof^ + ^30 Jaf + ^30 t/aj''^^') . (2. 1) 

^^iijig/e NLO subtraction term for 4-parton configurations, daf'^^ and da^'^^ are generic 

NNLO subtraction terms, which can be further decomposed into 

d4'''^ = da[""P + daP'"'''"' - dcc^'""'" - daf + ^/af (2.2) 

In a hybrid scheme of subtraction and slicing the subtraction terms have to satisfy weaker condi- 
tions as compared to a strict subtraction scheme. It is just required that 



(a) the explicit poles in the dimensional regularisation parameter e in the second line of eq. (^Jj) 
cancel after integration over unresolved phase spaces for each point of the resolved phase 
space. 



(b) the phase space singularities in the first and in the second line of eq. (2. 1 ) cancel after azimuthal 
averaging has been performed. 

Point (b) allows the determination of the subtraction terms from spin-averaged matrix elements. 
The subtraction terms can be found in [^-^]. The subtraction term daf'^^ without da'^^'^' 
would approximate all singularities except a soft single unresolved singularity. The subtraction 
term da"'^' takes care of this last piece |2^. The azimuthal average is not performed in the 
Monte Carlo integration. Instead a slicing parameter T] is introduced to regulate the phase space 
singularities related to spin-dependent terms. It is important to note that there are no numerically 
large contributions proportional to a power of Inrj which cancel between the 5-, 4- or 3-parton 
contributions. Each contribution itself is independent of Tj in the limit Tj — 0. 



3. Monte Carlo integration 

The integration over the phase space is performed numerically with Monte Carlo techniques. 
Efficiency of the Monte Carlo integration is an important issue, especially for the first moments 
of the event shape observables. Some of these moments receive sizable contributions from the 
close-to-two-jet region. In the 5-parton configuration this corresponds to (almost) three unresolved 
partons. The generation of the phase space is done sequentially, starting from a 2-parton config- 
uration. In each step an additional particle is inserted In going from n partons to n + \ 
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partons, the «+ l-paiton phase space is partitioned into different channels. Within one channel, the 
phase space is generated iteratively according to 

d^n+\ = d^nd^unresolved i,j,k (3.1) 

The indices /, j and k indicate that the new particle j is inserted between the hard radiators / and k. 
For each channel we require that the product of invariants SijSjk is the smallest among all considered 
channels. For the unresolved phase space measure we have 

1 1 2n 

d$unresolvediJ,k = J djCi J dX2 j dip @{\ - Xi - X2) (3.2) 



We are not interested in generating invariants smaller than (tj^), these configurations will be re- 
jected by the slicing procedure. Instead we are interested in generating invariants with values larger 
than (t]^') with a distribution which mimics the one of a typical matrix element. We therefore gener- 
ate the + l)-parton configuration from the ?i-parton configuration by using three random numbers 
Ml, M2, M3 uniformly distributed in [0, 1] and by setting 

^i = '^ps-' ^2 = 'nps (p = 2nu^. (3.3) 

The phase space parameter rjps is an adjustable parameter of the order of the slicing parameter rj. 
The invariants are defined as 

Sij=XiSijk, Sjk=X2Sijk, Sik = {l-Xi-X2)Sijk. (3.4) 

From these invariants and the value of (p we can reconstruct the four-momenta of the {n + l)-parton 



configuration [24|. The additional phase space weight due to the insertion of the (n + l)-th particle 
is 

Note that the phase space weight compensates the typical eikonal factor Sijk/{sjjSju) of a single 
emission. As mentioned above, the full phase space is constructed iteratively from these single 
emissions. 

4. Numerical results 

Fig. |l] shows the results for the Durham three jet rate and the thrust distribution at the LEP 
I centre-of-mass energy = mz with as{mz) = 0.118. The LO, NLO and NNLO predictions 
are shown together with the experimental measured values from the Aleph experiment [^]. The 
bands give the range for the theoretical prediction obtained from varying the renormalisation scale 
from jj. = Q/2 to jJ. = 2Q. Note that the theory predictions in these plots are the pure perturbative 
predictions. Power corrections or soft gluon resummation effects are not included in these results. 

In a recent calculation the logarithmic terms of the NNLO coefficient of the thrust distribution 



have been calculated based on soft-collinear effective theory [ ]26| ] 



—L = -\a^\n^T + a4\n^z + a3lrv'T + a2ln^'T + ai\nT + ao + <ff(T)], z = l-T. (4.1) 
dz X 
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Durham three -jet rate 



Thrust 




Figure 1: The scale variation of the Durham three jet rate and the thrust distribution at = mz with 
OJi (niz) = 0.118. The bands give the range for the theoretical prediction obtained from varying the renor- 
mahsation scale from jti = mz/2 to fi = 2mz. 
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Figure 2: A comparison of the NNLO coefficient of the thrust distribution as obtained from the numerical 
program with the logarithmic terms obtained from SCET. 



The values of the a/s are for Nf = 5 

a5 = -18.96, a4 = -207.4, ^3 = -122.3, = 1488.3, ai = -822.3, ao = -683.4. 

The logarithmic terms give a good description of the thrust distribution in the close-to-two jet 
region. They are not expected to give an accurate result in the hard region. Fig. ^ shows the com- 
parison of the NNLO coefficient of the thrust distribution as obtained from the numerical program 
with eq. (p~T|). In the left plot of fig. ^ the x-axis shows (1 — T) on a linear scale. This corresponds 
to the hard region, where the NNLO result from the numerical program is expected to give the cor- 
rect answer. The middle plot of fig. § shows (1 — T) on a logarithmic scale around {I — T) ^ 0.1. 
This corresponds to the peak region or the overlap region, where the perturbative NNLO result and 
the one obtained from SCET agree. The right plot of fig. ^ shows (1 — T) on a logarithmic scale 
around (l — T)^ 0.001. This corresponds to the extreme two-jet region, in which the logarithmic 
terms are dominant. In this region the results from the numerical program show a dependence on 
the slicing parameter. The numerical results for T] = 10~^, T] = 10~^ and T] = 10~^ are plotted. For 
smaller values of 7] the SCET result is approached. 
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